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■ Abstract. The traditional Hamiltonian structure of the equations governing 
conservative Rayleigh-Benard convection (RBC) is singular, i.e. it's Poisson bracket 
possesses nontrivial Casimir functionals. We show that a special form of one of these 

■ Casimirs can be used to extend the bilinear Poisson bracket to a trilinear generalised 
| Nambu bracket. It is further shown that the equations governing dissipative RBC can 

c/2 ' be written as the superposition of the conservative Nambu bracket with a dissipative 

^ , symmetric bracket. This leads to a Nambu-metriplectic system, which completes the 

geometrical picture of RBC. 
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The noncanonical Hamiltonian form of the hydro-thermodynamical equations in 
Eulerian variables is typically singular. This gives rise to the existence of a special 
class of conserved quantity, the Casimir functionals. This singularity is a consequence 
of the reduction that takes place if one changes the coordinates from the (canonical) 
Lagrangian coordinates to the (noncanonical) Eulerian coordinates by means of the 
particle-relabeling symmetry (e.g. Thiffeault and Morrison 2000). 

In the noncanonical Hamiltonian form the existence of additional conserved 
quantities is still hidden and hence it is natural to seek for a formulation which lets 
enter them the description in a similar way as the Hamiltonian does. For the great 
majority of hydrodynamical systems this is possible (Nevir and Blender 1993, Nevir 
1998) and leads to a form of description that formally resembles the structure of Nambu 
mechanics, which was first introduced by Nambu (1973) for discrete systems. 

We start with the classical definition of a discrete Nambu system (Takhtajan 1994): 

Definition: Let M be a smooth manifold and C°°(M) the algebra of infinitely 
differentiable real valued functions defined on M. Then M is called a Nambu-Poisson 
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manifold of order n if there exists a map {,...,} : [C°°(M)]® n — > C°°(M) called the 
Nambu bracket that satisfies 

(i) {fi, f2, ■ ■ ■ , fn} = (-iy ((7) {fa(i)J*(2), ■ fa(n)} where a denotes an element of the 
symmetric group of n elements, with e(a) being it's parity. 

(h) {/1/2, fz, ■ ■ ■ , fn+l} — /l{/2, f$,---, fn+l} + {fl, fz, ■ ■ ■ , fn+l}f2 

( m ) {{fl, • • • , fn-1, fn}, fn+l, ■ ■ ■ , /2n-l} + {/n, {/l • • • , fn-1, fn+l}, fn+2, ■ ■ ■ , /2n-l} + 
1" {fn, ■ ■ ■ , fln-2, {fl, ■ ■ ■ , fn-1, /2n-l}} = {fl, ■ ■ ■ , fn-1, {fn, ■ ■ ■ , /2n-l}} 

Note that the first property defines the skew-symmetry of the Nambu bracket, the second 
property is the Leibniz rule and the last one is a generalisation of the Jacobi identity, 
known as the Fundamental identity (FI) or Takhtajan identity. 

The Nambu bracket defines the kinematic part of a Nambu system, which is 
supplemented by n — 1 function Hi, ... , H n _i, such that the evolution of a real- valued 
function on M is given by 

-77 = {f, Hi, . . . , i? n -i}- 

Due to the antisymmetry of the Nambu bracket, it follows that Hi,...,H n _i are 
conserved by the Nambu flow. 

The generalisation to field equations was done by Nevir and Blender (1993) 
and Nevir (1998) by starting with the noncanonical Hamiltonian formulation of the 
respective model equations. They extended them by means of using one of their Casimirs 
as additional conserved quantity. Hence, their continuous Nambu formulation also only 
uses trilinear bracket structures. This kind of generalisation is the one we aim to use 
for the equations governing RBC. 

Let us now turn to dissipative systems with a conservative Hamiltonian core. Although 
the notion of such so called metriplectic systems is not unique (see Guha 2007 for a 
review) the constituent parts of them are the antisymmetric Poisson and a symmetric 
(or gradient) structure, that accounts for dissipation. Adding both pieces together then 
describes the dynamics of the whole dissipative system. The metriplectic bracket hence 
reads 

(if, 9)) :={f,g} + (f,g) 

where 

{/, g} = —{g,f} antisymmetric (Poisson) bracket 

(/, g) = (g,f) symmetric (gradient) bracket 

holds. 

In this paper we show that the equations of RBC can be cast in Nambu-metriplectic 
form, that is, the conservative part possesses a Nambu and the dissipative part a 
symmetric bracket structure. 
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2. The Nambu structure of conservative RBC 



The equations of two-dimensional RBC in case of an incompressible fluid using the 
Boussinesq approximation in nondimensional form read (e.g. Thiffeault and Horton 
1996): 



As usual, ip is the stream function generating two-dimensional nondivergent flow 
in the x-z-plane, ( = V 2 ip is the vorticity, T is the temperature departure of a 
linear conduction profile, v and k are kinematic viscosity and thermal conductivity, 
respectively. [a,b] := da/dxdb/dz — da/dz db/dx denotes the Jacobian operator. 

In the following, we will assume that the domain of interest is periodic in the re- 
direction. In the vertical, we either also assume periodicity or stress-free boundary 
conditions (Saltzman, 1962), i.e. 



^ = o, C 



^ 

' dx 







at 



0, 



In either case, since in RBC a constant temperature difference is maintained externally 
between the top and bottom of the fluid, the appropriate boundary condition for the 
temperature deviation T is 



T = 



at 



z — 0, z — z too . 

Determining the classical continuous Nambu form is possible only in the case of vanishing 
dissipation, i.e. in case v — k — 0. Then, since the remaining terms on the right hand 
side can be written as dT/dx = [T,z] and dip/dx = [if),z], respectively, both equations 
may be arranged as 



f t +[^,Q + [z,T-z] 



0, 



dT 



0. 



(2) 



This set of equations is Hamiltonian upon using 



« = / o (I W )'-T,)d/ 



as Hamiltonian functional and 
5 J 7 5H 



C 



+ (T-z) 



df 



(3) 



as a Poisson bracket. An analogue bracket arises also in magnetohydrodynamics (MHD) 
(Morrison and Hazeltine 1984). Note, however, that this Poisson bracket is singular, i.e. 
it possesses nonvanishing Casimir functionals. They are 

Ci = / g(T - z) df, C 2 = [ (h(T - z) df, 

Jn Jn 

where g, h are arbitrary functions of T — z. For a physical intepretation of the analogue 
Casimirs in MHD, see Thiffeault and Morrison (1998). In conservative RBC the first 
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class of Casimirs C\ physically describes the preservation of T — ^-contours. In fact, this 
class of Casimirs generally arises in nondivergent and inviscid fluid models. In turn, the 
second class C2 incorporates the Kelvin's circulation theorem, which in the case of RBC 
requires conservation of the vorticity ( on closed contours of T — z. 

To determine the conservative Nambu form, we are only interested in the special 
form of C2 with h = T — z: 



C 



C(T - z) df. 



It allows us to extend the Poisson bracket formulation of the governing equations to 

= {(,C,H} 

-{T,c,n}, 

■} is defined for arbitrary functionals J-%, T<i and JF 3 by the 



arrive at the their Nambu form: 




d( 


~5C 5H 




~5C 5H 


dt 


_5T , J(_ 






dT 


~8C 5H~ 




~dt~~ 







where the bracket {• 
equation 

\T 1, JF 2 , T 3} := - 



5F 2 5^3 



+ 



8j\ 



5F 2 5-^3 



5T' 5( 



+ 



8F1 



8T2 <5JS 



5( ' ST 



df- (4) 



This bracket is easily seen to be totally antisymmetric in case we assume periodic 
boundary conditions in both directions. Another possibility that guarantees 
antisymmetry of (jl]) is provided by periodicity in x-direction and free boundaries in 
the vertical. For this choice, however, it is necessary to fix JF 3 = Ti in the Nambu 
bracket. This is explicitly shown in the appendix. Assuring antisymmetry is also 
possible by only considering those class of functionals that sufficiently rapidly go to 
zero towards the boundaries. All these choices guarantee that surface terms emerging 
from an integration by parts vanish. 

The FI is proved by noting that the above Nambu bracket is simply the continuous 
analogue of the heavy top Nambu bracket. Hence, this bracket may indeed serve as a 
good Nambu bracket. 

Note that C is indefinite with respect to sign. In this respect it is akin to the helicity, 
which is a Casimir for the three dimensional incompressible Euler equations (Nevir and 
Blender, 1993). To our knowledge, the latter model and RBC are the only known ones 
that need indefinite Casimirs to allow for a Nambu representation. 



3. The Nambu-metriplectic structure of dissipative RBC 

Now turning to the dissipative equations (i.e. v 7^ 0, k 7^ 0), we aim to show that this 
model has Nambu-metriplectic form. Let us first note, that adding = [T — z,T — z] 
to the first equation in ([2]) leads to the equivalent system 

§ + [V, C] + [T, T - z] = 0, ^ + [V, T — z] = 0. 
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which is Hamiltonian upon using 

1 



as Poisson bracket and 



-(V^) 2 -Tz--(T-z) 2 



df = H-S 



as Hamiltonian functional. Note that S is the particular realisation of the first class of 
Casimir functionals C\ with g = 1/2(T — z) 2 . In Morrison (1986) functionals like Q are 
termed generalised free energy. 

This formulation enables us to cast the equations governing dissipative RBC in 
Nambu-metriplectic form. Indeed, in this case the equations can be written as 



d_C + 
dt 


'5C 


6G~ 


+ 


~5C 5G~ 


Jf 1 


sc. 




dT 

a + 


~5C 

Jc 


sg~ 

sc 




ST 



-z/V 



j_g 



and it is possible to introduce an indefinite symmetric bracket that governs dissipation: 

(5) 



The symmetry property of this bracket is assured if either periodic boundaries are 
assumed or only functionals that sufficiently rapidly decay to zero near the boundaries 
are considered. Adding together both brackets gives the entire dynamics of two- 
dimensional RBC: 



§ = {C,C,f?} + <£0> 



^ = {T,c,g} + (T,g). 



Note that in some sense the Nambu-metripectic formulation of (CQ) is geometrically the 
most complete, since representatives of both classes of Casimirs are needed to represent 
the whole dynamics. 

The existence of dissipation generally spoils the conservation properties of 
conservative system. This is also the case in RBC, since we have 



f = <C,0>*0. 



That is, the evolution of g is determined solely by g, whereas the evolution of C is 
determined both by C and g. 



4. Comments and outlook 

In this work we have shown that the conservative part of (pQ) can be written in Nambu 
form, while the full dissipative system possesses a Nambu-metriplectic form. For both 
representations, the explicit usage of Casimir functionals of ([3]) is crucial. The Casimir C 
allows to extend the bilinear Poisson bracket to a trilinear Nambu bracket. In turn, the 
Casimir S can be subtracted from the Hamiltonian 7i to give the modified Hamiltonian 
g. This doesn't alter the dynamics, since Casimirs are trivial conserved quantities 
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and Hamiltonians are only determined up to Casmir functionals. But introducing Q is 
essential to allow for the necessary symmetry property of the dissipative bracket (jSJ). 

The Poisson bracket ([3]) is an example of a Lie-Poisson system. Such systems are 
built from an underlying Lie algebra structure, typically owing to a reduction from a 
set of canonical to a set of noncanonical variables (e.g. Thiffeault and Morrison 2000). 
In case of RBC it is the semi-direct extension of the algebra associated to the group 
of volume-preserving diffeomorphisms on some domain Q with the vector space of real- 
valued functions on Q (Thiffeault and Morrison 2000). As was shown in our work, the 
Nambu structure is compatible with this type of algebra extension. In particular, the 
Nambu bracket (jlj) allows to put ([3]) in a more symmetric form. 

Studying algebras and extensions thereof in conjunction with Nambu structures 
offers a way to classify Nambu systems. Such work is currently in progress and will be 
published elsewhere. 
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Appendix 



We aim to explicitly show here that assuming periodic boundaries in x-direction and 
free boundaries in the vertical guarantees the total antisymmetry of (jlj) provided we fix 
the Hamiltonian 7i in the bracket. 

We only need to prove antisymmetry with respect to the first two functionals of the 
bracket (jlj). The antisymmetry in the other pair of arguments is obvious. An integration 
by parts gives 

f 6C \5JF_ 8H 

6C_ \5JF_ sn 
, n 5T [JC'JC 
SC \5JF_ 5H 

were we have already taken into account periodicity in the x-direction. Due to the 
imposed boundary conditions 



a *C 
8JF_ 

a *C 



~6C 5H 
6C_ 5H 

6C m_ 

5C' Jf 



df 
df 

df- 



df- 
df- 
df- 



/6^6Cd_6H 
L \6T J( dx Jc 

l \KVTdx~JCj 

L \6( Kdx'JT) 



dx 



dx 



dx 







dip 


dx 


If) 


z dx z 



0. 







dz 


dx 


m 


z dx 







all the second terms on the right hand side vanish. Similar considerations also hold for 
fixing the modified Hamiltonian Q. 
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